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jj;^ I Abstract 

^ . We sketch a group-theoretical framework, based on the Heisenberg-Weyl group, en- 

! compassing both quantum and classical statistical descriptions of unconstrained, non- 

' relativistic mechanical systems. We redefine in group-theoretical terms a kinematical 

^ ■ arena and a space of statistical states of a system, achieving a unified quantum-classical 

! language and an elegant version of the quantum-to-classical transition. We briefly discuss 

I the structure of observables and dynamics within our framework. 

(N ; 

^ 1 Introduction. 

00 ■ ^ ^ 

O ■ Since the seminal works of Weyl [I] and Wigner [2], the fundamental role of group theory 

^ . in quantum mechanics has become an established fact. The aim of the present work is to 

\^ [ remark that it is possible to redefine, using solely group-theoretical notions, the mathe- 

O ' matical representations of kinematical arena and state-space in non-relativistic quantum 

^ ■ mechanics, such that i) a unified language for quantum and classical statistical descriptions 

Oh! is obtained; ii) there is a natural transition mechanism, leading from a more generic quan- 

I tum state-space to a classical one. The latter is a mathematically precise formulation of 

■ the classical limit of quantum theory at the level of statistical descriptions. It constitutes 
^ , the main result of our paper. Since the literature on the subject is enormous, let us state 
^1 it clearly that our work is neither meant to be an overview of the quantization methods, 

■ nor even an overview of various realisations of classical limits of quantum mechanics. For 
^ , that see, e.g., Refs. [3l|4| and the references therein. 

^ I In classical physics, statistical description of mechanical systems is given in terms of 

- - ' probability measures on a phase space T. On the other hand, in quantum mechanics, 

statistical properties are encoded into density matrices g acting on a Hilbert space TC. 
A natural question arises how to connect these two such seemingly different formalisms 
through a sort of a quantum-to-classical transition. Obviously, such a transition must 
exists, as indirectly proven by numerous examples and more or less heuristical arguments 
within all possible approaches to quantum theory. The question is rather how to write it in 
a clear, precise manner in a hope to shed some light on how the classical world appears. One 
strategy, which we adopt in this work, is first to try to find a common theoretical framework 
for both quantum and classical theories and then search for a transition mechanism within 
this framework. The standard realization of such approach is to use Wigner [5] or Moyal 
functions |6j, associating with every density matrix g a phase-space pseudo-probability 
distribution Wg and a corresponding pseudo-characteristic function Xg respectively (e.g. 
Ref. |7] contains a modern exposition). However, it is a well known fact [HIT] that both 
Wg and Xg fail to satisfy the positivity conditions that possess their classical analogs. 
Moreover, there seems to be no direct relation between positivity of Wigner or Moyal 
functions and "classical behaviour" of density matrices. Whether this is a drawback or not 
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in the context of providing a unified framework for both quantum and classical statistics 
is perhaps a matter of taste, but let us assume that it is. 

The aim of our work is to point out that such a unified language with a clearly visible 
quantum-to-classical transition mechanism is achieved if instead of working with den- 
sity matrices and probability measures one performs (generalized) Fourier transform and 
works with (generalized) characteristic functions. The structure, which emerges after such 
Fourier transform is universal — in both quantum and classical cases it consists of a certain 
group G, serving as a sort of "kinematical arena" (instead of a phase space F or a Hilbert 
space 7i) and a set of normalized, positive-definite functions (p, representing statistical 
states. Recall [8] that a complex function i;^ on a group G with a Haar measure dg is called 
positive-definite if it is bounded, continuous, and satisfies: 

J [ dgdhj{g)4>{9-^h)f{h) > for any / G L^G); (1) 

by normalization we mean here that: 

0(e) = 1. (2) 

In the case when G = M", Eqs. ([T| and jj]) represent familiar properties of a characteristic 
function of a classical probability distribution fj. on |^. As we shall show in the sequel, 
when G is non-Abelian, cp can still be viewed as a characteristic function, which we call 
"non-commutative characteristic function", but of a quantum probability distribution g. 
And this is precisely the difference between classical and quantum statistics in the emerging 
formalism — in quantum case group G is necessarily non-commutative, whereas in classical 
case it is Abelian. 

In order to avoid the need for performing Fourier transform forth and back and thus 
achieve some, at least formal, conceptual simplicity, we propose to make one step further 
and postulate that a suitable group G together with a (sub)set of normalized, positive- 
definite functions (p should be taken as a basis of statistical description of mechanical 
systems, both classical and quantum. Thus, the change of concepts we propose to exam- 
ine is the following: i) as the kinematical arena of the statistics we consider a certain 
group G, called "kinematical group", together with its irreducible unitary representations; 
ii) as mathematical representatives of statistical states of the system we consider normal- 
ized, positive-definite functions on G. We show that such a unified formalism is indeed 
equivalent to the standard ones, upon a correct choice of G. 

The kinematical group of our approach should not be confused with a group of symme- 
try transformations of neither kinematics nor dynamics of the theory (for an alternative 
programme for quantum theory, where G is taken to be a group of dynamical symmetries 
see Ref. flQ]). Its role is rather to serve as a background for statistics (just like classi- 
cal phase space F or Hilbert space H), encoding statistical properties of the system in a 
proper way. For example, as we show in this work, in the case of unconstrained mechan- 
ical systems (particles) with n degrees of freedom the right kinematical group turns out 
to be the Heisenberg-Weyl group Hn [HIUlIl]. Intuitively, this can be understood in the 
following way: since in the standard formulation of quantum mechanics Hn is a realization 
of Heisenberg uncertainty principles, it turns out that it is this information that is enough 
to re-produce quantum statistics. Hence, we invert the usual role of the Heisenberg-Weyl 
group and instead of treating it as a mere consequence of the uncertainty principles, we 
propose to look at it as the source of the latter. In this sense our approach bears some sim- 
ilarities to Klein's Erlangen programme — as we shall see in the sequel kinematical group 
determines convex geometry of statistical states. However, let us stress again that, unlike 
in Erlangen programme, our kinematical group is not a group of symmetries. 
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When viewed from the perspective of the standard formahsm, the resulting approach 
closely resembles Weyl quantization [Ij. The difference is that instead of quantizing func- 
tions on the classical phase space F = R^", which can be viewed as an Abelian group, we 
quantize functions on the non- Abelian group Hn ~ M x F. This seemingly subtle difference 
produces, as we will show, some interesting changes. For example, the Hilbert space struc- 
ture of quantum theory does not have to be fully and independently postulated, but to 
some extent follows naturally from the consideration of irreducible representations of Hn, 
if the latter is treated as a fundamental entry of the formalism. Second, and perhaps more 
importantly, one clearly sees how the classical statistical state-space naturally emerges 
from the quantum one as Hn collapses to one of its Abelian subgroups. Thus, the essence 
of quantum-to-classical transition in our framework is the restoration of commutativity 
of the kinematical group. In the context of more modern versions of quantum theory, 
our formalism is closely related to the algebraic approach (see e.g. Ref. [12] for a deep 
exposition and Ref. [4] for the latest trends) and, in fact, can be viewed as a concrete, but 
rather non-standard, realization of the latter. 

Finally, let us mention that in physical literature, the group-theoretical formalism that 
we develop in the present paper was in fact initiated by Gu in Ref. [13]. Especially 
in the context of providing a more coherent, as compared to the standard Wigner and 
Moyal functions, way of describing both quantum and classical statistics. However, Gu 
did not fully perform the reformulation of the theory and concentrated mostly on practical 
problems, treating non-commutative characteristic functions rather as secondary objects 
with respect to the usual density matrices. Neither did he examine the representation of 
observables and the classical limit (on the level of kinematics) in the resulting formalism. 
In the present work we explicitly carry over the mentioned conceptual change and treat 
from the beginning non-commutative characteristic functions as primary objects of the 
theory, while density matrices or probability measures as secondary. 

The plan of the work is the following: in Section [2] we recall the basic properties of 
the Heisenberg-Weyl group. In Section [3] we develop the group-theoretical formalism and 
present our version of the quantum-to-classical transition. In Section [4] we sketch the 
group-theoretical representation of observables and briefly comment on the dynamics in 
our scheme. Then, in Section [5l we show with two physical examples how the classical 
transition in our language works in practice. It should be kept in mind however that the 
objective of this work is not to develop practical methods of calculating classical limits of 
density matrices, but rather to put quantum-to-classical transition in a mathematically 
rigorous form. The concluding remarks are gathered in Section [6l 



2 Heisenberg-Weyl group 

The basic object of our study will be the Heisenberg-Weyl group Hn, where n is the number 
of degrees of freedom of the considered mechanical system. Thus, at this moment we fix 
the kinematical group: G = Hn- The group Hn can be identified with a space M x M" x R", 
equipped with the following multiplication law: 

{s, ri, ^) ■ is', r?', ^') := (^s + s' + iu;[(r,, iv',^%r] + ri',^ + ^'), (3) 

where {s,r],^) are the coordinates and: 




In the sequel we will interchangeably denote group elements by g,h, . . . or by the corre- 
sponding coordinates. The Haar measure dg on Hn is just dsdL^^d^rj. 
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In what follows we will need irreducible, unitary (and strongly continuous) represen- 
tations of Hn- They are characterized by the Stone-von Neumann Theorem (see e.g. Refs 
[8j E]). Let us briefly recall their structure. There is a family of infinite-dimensional 
representations , M. 3 X 0: 



-iXs 



exp 



(5) 



(throughout the work the repeated indices are summed over), where the self-adjoint gen- 
erators q, p satisfy on the common domain: 

[qj,Pk] =iX6jk, (6) 

and a family of one-dimensional representations Tq p, labelled by points (q, p) G M^": 

Tlp{s,v,^):=e'^^^'i^-^^P^K (7) 



From experiment we know that the representation realized in Nature is with X = h. 



3 Group-theoretical approach 

In order to expose the quantum-to-classical transition mechanism, we adopt the strategy 
that first a suitable unified quantum-classical statistical framework should be developed. 
Traditionally, a search for such a framework has been interpreted as a search for quantum 
analogs of classical probability distributions. In case of mechanical systems one thus 
follows Moyal [6], and defines a phase-space characteristic function (also known as the 
Moyal function) corresponding to a given density matrix g by: 

Xeiv, ^) := tr(^Qexp ^{VjQj - CjPj) ) • (8) 

One passes then to its Fourier transform — the Wigner function [5l[6]: 

/d"ifd"r) 
^^l^e-i(v.i.~im)^^^rj,^), (9) 

in the hope to obtain an analog of a classical probability distribution. However, this 
attempt fails — as we have mentioned in the Introduction, Wigner function jO]) is generically 
non-positive on the classical phase space F = R^". Moreover, there seems to be no universal 
relation between positivity of Wg and "classical behaviour" of density matrices: there are 
density matrices showing what is generally accepted as "genuine quantum behaviour", and 
nevertheless possessing positive Wigner functions (for example, so called, squeezed states 

m) 

There has been developed some methods to get around the above difficulty. One of 
them is to replace the abstract definitions ^ and ([9|) by operational ones, i.e. involving 
prescribed interaction with an external reference particle. This allows one to construct a 
positive phase-space probability distribution (see e.g. Ref. [HI). Another way of producing 
a positive phase-space probability distribution is to use the Glauber-Sudarshan coherent 
states I a) [IHIIT], and assign to each density matrix g the Husimi function, also known 
as the Q-representation, {a\ga) (see e.g. Ref. [I6]). 

In contrast to the approaches mentioned above, we propose, following Ref. [13], an 
alternative way towards the unification of languages of quantum and classical statistics. 
Instead of searching for positive phase-space probability distributions for density matri- 
ces, let us rather change the object of our interest and look at the properly generalized 



4 



characteristic functions. Observe that the lack of positivity of Wg{q, p) is mathematically 
related to the fact that characteristic function ^ generically is not positive-definite (c.f. 
Eq. ([I|)) on R^" |T7|. Gu noted in Ref. [13] that if one replaces the standard definition 
^ by, in our opinion, more logical one: 



then such defined function (j) is positive-definite on Hn and is also normalized: (/)(e) = 
1, where e = (0,0,0) is the neutral element. Hence, as we have mentioned earlier, it 
possesses all the features of a classical characteristic function of a probability distribution. 
However, unlike the latter, (pg is defined on the non-Abelian group H^, rather than on E?^. 
This justifies the terms "kinematical group" for Hn and "non-commutative characteristic 
function" for (/)g. Note that one cannot straightforwardly apply the usual Fourier transform 
to (pg, as it was done in Eq. ([9]). Rather Eq. is a non-commutative Fourier transform 
of the density matrix g. 

Motivated by the above observation, we propose to examine the following alternative 



construction of quantum statistics of an unconstrained mechanical systenu: 

• treat the group Hn as the basic entry of the formalism, which sets up the kinematical 
arena; 

• take as statistical states of the system (abstract at this moment) normalized, positive- 
definite functions (j) on Hn] the set of such functions will be denoted by Vi{Hn) and 
it is a convex subset of the set of all continuous, bounded functions on Hn- 

Heuristically, the appearance of Hn rather than F may be explained as follows: Hn = 
M X M^" ~ M X F, so we may view (to some extent) Hn as an extension of the classical 
phase space. The additional degree of freedom, labelled by s, per point of F can be then 
attributed to quantum-mechanical phase, which is supported by the multiplication law dH]) 
and the form of generic representations jS]). This phase degree of freedom is generally non- 
compact (it is diffeomorphic to M), but from Eq. ^ we see that once we work in a fixed 
representation of H^ which we will do in what follows, it effectively becomes U{1). This 
makes our formalism loosely resemble Kaluza-Klein theory [19] — an approach to unification 
of Maxwell and general relativity theories. There one extends space-time through adding 
(in a local way) a phase degree of freedom at each space-time point. By introducing 
a suitable parallel transport on this 5-dimensional space and postulating an analog of 
Einstein equations one then recovers coupled gravitational and Maxwell fields. Of course 
the analogy is only distant, as we are not introducing a parallel transport on our extended 
phase space, but rather the group structure (cf. Eq. jS])). The motivation for the choice 
of states is more straightforward — in both in classical and quantum cases characteristic 
functions possess the same features, provided they are properly defined through Eq. ([TO]) . 

However, the whole Vi{Hn) turns out to be too large. This happens because a generic 
(j) G Vi{Hn) contains contributions from all possible representations of Hn, while we know 
that only one of them is realized in Nature. To identify the set of physically relevant 
states within Vi{Hn) and recover the standard density matrix formalism, we use Gelfand- 
Naimark-Segal (GNS) construction (see e.g. Ref. [8]). Although this is the fundamental 
tool in algebraic approaches to quantum theory, note that here we are using it in a different 
manner. In particular, we are not starting from a C*-algebra of observables, but rather 
from the kinematical group [20]. Using the GNS construction we can uniquely (up to a 
unitary transformation) assign to each abstract state 4> G Vi{Hn) a triple (W^, vr,^, f,/,), 
where tt^ is a representation of Hn acting in a Hilbert space 7i^, is a normalized cyclic 

For a similar approach, based on covariance systems see Ref. [18] 



^g{g) ■.= tr[QT\g)]=e-'^Xe{ri.$.), 



(10) 
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vector, and: 



(t>{9) = {v^\T^^{g)v^) . (11) 

The representation tt,^ is generically reducible. It is irreducible if and only if (j) is an 
extreme point of the state-space Vi{Hn). Such states will be called "pure". Motivated by 
the commutation relations ([6|), we call physical states those abstract states, for which tt^ 
is a countable multiple of T^: iTfj, = • T^, TC^ = 0^ Tin., since then: 

m = {v^\®T%H) = J2(''^\T^(9)v^) = tr[0^T%)], (12) 

i i 

where fj's are the components of in each copy of Hn, and: 

e-fli := X]^*lTrT)'vTrTl I ' Pi-=\\'"i\\^^ = \\v<i>\\'^ = 'i- ■ (13) 

Hence, to each physical state 4> we may assign a positive trace-class operator 'q^ in 71%, 
representing (/>. Prom Eqs. ^ and (fT2]) we infer that physical states are of a special form: 

0(^,^7,0 = e-'^x(r?,0, (14) 

where x is simply the standard characteristic function ([8|) of ^(^. 

The converse also holds, i.e. each abstract state 4> of the form ([14| is physical and we 
can uniquely assign to it a density matrix This has been in fact proven in Ref. [13j. 
The operator is defined as follows: 



• I ^"^''"\mT"(,)\ (15) 



„ (2„f! (2ifi) 

(provided the above integral exists in the sense of matrix elements, for which it is sufficient 
that X G L1(]R2«)). Since tr(^exp ^ir^^q^ - ^jpj) ) = (27r?i)" 5"(^)(5"(77), we have that: 

<P,,{g)=tr[Q^T^{g)] =<t>{g), (16) 

and hence formula (fT5]l can be viewed as the inverse, with respect to the definition p^ . 
non-commutative Fourier transform. Prom Eq. (|16p and the uniqueness of the GNS 
construction, it then follows that g^f, is the same (up to unitary rotation) as the density 
matrix from Eq. (fT2]) . The representation (fTTjl is recovered by spectrally decomposing 
Qfj, and then going back from Eq. (fT2]l to Eq. (fTTl) . Moreover, if we look from the standard 
formalism point of view, then we also have g^^ = g (since matrix elements of satisfy 
orthonormality relations like matrix elements of an irreducible representation of a compact 
group; see Ref. |13|). 

Hence, physically relevant states are faithfully represented by functions (j) G Vi{Hn) of 
the form (fHlln . Please note that restriction to the above form does not break any group 
symmetry, since the kinematical group Hn is not supposed to act on the abstract state- 
space Vi{Hn) as a group of symmetries. This specific form only fixes the representation 
appearing in the GNS construction to the physically relevant one with X = h. Prom now 
on we will assume that we work only with the physical states. We have thus achieved the 
desired reformulation of the standard theory. Density matrices are now secondary objects. 



^In fact, for establishing this correspondence we could have used only the formulas (fTO|) and (fT5|) . but the 
GNS construction is more general — it can be carried out on an arbitrary locally compact group. 
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constructed from the physical states and the proper representation of the kinematical 
group through Eq. (fTSj) (or the GNS construction). We stress that we are deahng here 
with quantum statistics only, as the notion of the linear superposition seems not to be 
easily visible in the group-theoretical language. Nevertheless the linear structure pertains 
in our formalism through the GNS construction. It is then an interesting question weather 
some modification of the group-theoretical framework can be used for construction of a 
non- linear generalization of quantum mechanics. 

The main benefit of the presented reformulation lies, in our eyes, in that it provides a 
natural transition to the regime h = 0, in which one recovers classical statistics. Indeed, if 
we accept that what is experimentally available are density matrices (for example through 
the state tomography technique), then we are led to study the irreducible representations 
of the kinematical group in the case ^ = 0. As can be seen from Eq. jT]), the irreducible 
representations of Hn become in this regime effectively the irreducible representations of 
the Abelian factor-group Hn/{{s,0,0);s G M} = M^", parametrized by (r/,^). Thus the 
phase degree collapses and the kinematical group turns into the classical phase space. 
Hence, the states that we are naturally led to consider are now functions from Pi(M^"). 
We do not have to worry about fixing the right, physical representation, like in Eq. (fT4l) . 
as it is already fixed by setting h = - these are the representations ([7]). The crucial point 
is that Bochner's Theorem |8j states that the functions from Vi (R^") are in an one-to-one 
correspondence with Borel probability measures on the dual group M^", isomorphic to 
]^2n jgj rpj^g duality (•,•) is provided by the representation T^^ itself: ((q, p), (j?, ^)) := 

^q,p(^'^), (q>P) e M2n. Thus, in the classical regime our states, i.e. functions 4> G 

Vi{M?^), can be uniquely identified with Borel probability measures fj,^ on M^" ^ M^" and 
the latter space plays the role of the classical phase space of the system (at least in the 
context of statistical description). As a result, we recover classical statistical description 
of the system. 

In the case that (/> G Pi(M2«) is also in L1(M2"), or if we allow for distributions, we 
can explicitly recover ^u^ through the analog of the integral (fTSl) . which now becomes the 
usual Fourier transform (we only manually adjust the constant multiplying the measure): 

d/i<^ = <A(q,p) d"qd"p, (17) 
^(q>P) - / ^^(V,mlM = I ^^(^-Oe-^"^--^^-^), (18) 

and (/) is a classical probability density in the phase-space. 

When applied to composite systems, the above reasoning reveals one interesting aspect 
of quantum entanglement (see e.g. Ref. [21] for an introduction into the subject). Namely, 
when passing to the classical regime, the underlying kinematical group becomes Abelian 
and the corresponding classical states loose the ability to get entangled, since probability 
measures on Cartesian products can always be represented as suitable limits of convex 
combinations of product measures. Hence, the very existence of quantum entanglement 
may be linked, within our formalism, to the non-Abelian character of the kinematical 
group. More detailed study of the connection between group-theoretical methods and 
entanglement is discussed elsewhere |22| . 

4 Remarks on observables and dynamics 

In this Section we briefly remark on the representation of observables and dynamics in 
our group-theoretical language. We will not be very detailed and mathematically strict 
here, but rather present a general outline. The easiest observables to deal with are those 
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represented by trace-class operators in the standard language. In our reformulation they 
are given by complex continuous functions F from L^{Sn), Sn '■= [0, 27r//i] x M^" C Hn, 
satisfying: 



F{g-')=F{g). (19) 
The mean value of F in a state (p is defined as: 

(F)^:= / dg^{g)F{g), (20) 

Sn 

where we have rescaled dg so that it is now equal to dg = [(27r)^(27r/i)"~^] ^dsd'^^d"/;. 
The integral ([20]) is well defined due to the boundedness of cj). To establish the connection 
with the standard representation of an observable (for the connection to the algebraic 
approach see remark [20J), note that to each such F we can assign a hermitian operator 
Ap by a formula analogous to Eq. (fTSjl : 

Ap := f dgF{g)T\g). (21) 

J Sn 

The above integral exists, in the sense of matrix elements, as F € L}{Sn)- On the other 
hand, to each trace-class observable A we can assign a continuous function Fa by an analog 
ofEq. (fTOjl: 

FA{g):=iv[AT\g)\ (22) 

Using the same arguments as in the case of density matrices, one can easily show that 
Fap = F and Apj^ = A, thus establishing the correspondence between trace-class observ- 
ables and functions satisfying Eq. (fT9]l . 

In the case of observables not representable by trace-class operators, one has to allow 
for distributions. We will not investigate here which exactly distribution space one needs 
to consider in order to cover all relevant observables, but only write down the distributions 
Fq. and Fp. representing the generators qj and pj-. 

2 

Fq, {s, rj, ^) = ih {27rhrew [i^ - ^] '5"(^) • • • 5,,5(r/, ) . . . 5(r?„,) (23) 

r eh 

Fp, is, V, ^) = -i^ (27rn)«exp [ihs - 5^{ri) 6(^1) . . . %5(e, ) . . . di^n). (24) 

Higher order polynomials in qj and pj are proportional to higher order derivatives of 
Dirac's delta. 

As we have mentioned in the Introduction, the above procedure of recovering standard 
formalism from the group-theoretical one closely resembles Weyl quantization. Very briefly, 
Weyl proposed to assign with each Fourier-representable function F on the classical phase 
space r 

F(q,p) = j d''id'''nF{r],i)e^^''^'i^-^^P^^ (25) 
an operator 

j^Weyi _ y"d"^d"T7^(^,0e^^''^*"^^^^^ (26) 



A comparison of Eq. (|26]l with Eqs. ([21]) and (11511 shows that we basically substitute F 
with Hn in the original Weyl formalism. As we have shown in the preceding Section, this 
substitution leads to some concrete benefits. 



8 



Let us move to the representation of dynamics in our language. The dynamical law 
takes a form of a differential equation imposed on a path t (pt, which should be equivalent 
to the von Neumann equation for the corresponding density matrix Qt. In fact, this 
equation was derived and analyzed in Ref. [13] and we merely quote it here: 



H[- ihdr, + ihd^ + ^r?) -h{- ihdr, - ihd^ -^v] <Pt, (27) 



where H{q,p) is the Hamiltonian, which for simplicity we assume to be of the form 
H{q,p) = T{p) + V{q) with T and V analytical. Note that due to the property ([HI, Eq. 
([271) is, modulo the phase factor, just the quantum Liouville equation [7], but imposed on 
the characteristic function xiVi^) rather than on the Wigner function VK(q, p). 

In the classical limit, as we have argued before, the kinematical group effectively col- 
lapses to the classical phase space F = M^" and observables become functions (or distri- 
butions) on r. The condition (fTOl) now reads: 



F{-r^,-i) = F{'n,i), (28) 
and from Eqs. ([2T[) and ([7} we obtain that: 

-^S_^F{rj,^)e~'(w.-^^P.). (29) 



Thus, observables correspond now to real (because of the condition ([28[l ) functions on the 
phase-space M^". Using Eqs. ([TTl) . ([TSl) . and ([29]) the state average ([20]l becomes simply 
the average of AF{q,p) with respect to the measure defined by the state in question 

(F)^ = I ^^^F{V,i)Hv,^) = I df,4q,p)AF{q,p). (30) 



To complete the picture, note that the classical limit of the dynamical law ([27[) for 
H{q,p) = p^ jlra + T^(q) reproduces (upon rescaling of r/,^ - see the next Section) the 
classical Liouville equation | ll3j . 



5 Examples of classical limits 

Here we briefly show with two physical examples how the procedure of moving to classical 
regime = works in practice. However, let us stress again that it is not the goal of this 
work to develop another tool for studying classical limits of quantum states, but rather 
to examine how the formalism of non-commutative characteristic functions leads to the 
more coherent quantum-classical language and the natural description of classical limit 
of quantum statistics. Since concrete examples of physically interesting states has been 
available to us in terms of density matrices q anyway, we have to start from them. From 
this perspective, our approach obviously brings nothing new to the standard methods of 
Wigner and Moyal functions, as can it be seen from Eq. ([T4[l . Thus, we repeat that the 
main goal of the present work is conceptual rather than practical. 

The prescription for taking classical limits is rather simple: use the basic formula ([TO]) 
to calculate <f)Q for a given matrix q [13]. Next, check if there exists, in the distributive 
sense, a limit lim;j_^o (point limits are too restrictive). If a state ^ is to possess a classical 
limit at all, we naturally expect that (lim^j^o 4>) ^ 'Pi(M^'^), or in other words (lim^^o 
should be a classical characteristic function. If that is the case, we can use the prescription 
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(|17H18|) to retrieve the corresponding probability measure. If not, i.e. lim^j^o does not 
exists, or is not in Pi(M^"), then the state in question does not possess the classical limit. 

To illustrate the procedure, let us first consider convex mixtures of coherent states, i.e. 
density operators g with a positive Glauber-Sudarshan P-representation [l5j: 



d//(a, a)|a)(a| 



(31) 



where a := -^^{q + ip) and // is a probability measure on the classical phase-space M^. 
We stress that we consider coherent states here purely kinematicaly, without any explicit 
or implicit relation to the dynamics. They are defined as the isotropic states minimizing 
the Heisenberg uncertainty relations, arising from Eq. ([6|) and their particular importance 
for quadratic Hamiltonians does not concern us here. 
Substituting Eq. pTl) into Eq. (fTOl) we obtain: 



-ihs 



exp 



1 

4h 



(32) 



At a first glance, the last term in the integrand in Eq. (j32|l does not seem to possess 
any meaningful distributive limit when ^ — > 0. That would be quite counterintuitive, as 
the matrices of the form (fSTI) show a classical-like behaviour: for example the averages 
of normally ordered observables are equal to the phase-space averages with respect to 
However, note that the parameters rj, ^ are just some arbitrary coordinates on the 
Heisenberg- Weyl group and we are free to re-scale them. Actually, the specific form of the 
operator in Eq. jS]) was motivated by the physical dimensional analysis (the argument 
of the exponential function should be physically dimensionless) and in order to recover the 
group multiplication law ([3|) one has to re-scale rj, ^ by h. From another point of view, this 
rescaling is suggested by the Wigner function which can be rewritten as follows: 



/d"^d"r) 

If we accept the above arguments, we are led to consider: 



h—^0 



(33) 



(34) 



instead of lim;j_^o </'(s, as the proper classical limit (compare to the methods of Ref. 
[TH]). Then from Eq. §^ we obtain that: 



(j)g{s,hri,n£) ji{r],C) ■= / d/x(g,p)e- 



(35) 



By Bochner's Theorem (or by an easy direct inspection) jj G Vi{M?) and obviously the 
corresponding probability measure is just the measure /i itself. 

As the next example let us consider the Fock states |m): g = \m){m\. Again, we 
consider them just as kinematical examples. From Eq. (fTOl) we obtain that: 



0e(s,7?,O 



1 

2h 



(36) 



where Lm{x) = ^^q(— l)^2;^/[(m — k)\{kl)'^] is the m-th order Laguerre polynomial. 
Just for illustration's sake, we will considered here a rather uninteresting limit ^ — > of the 
fixed Fock state m = const. Obviously, this limit does not have much physical sense, but 
from a purely formal point of view the vectors \m) are legitimate states in the kinematical 
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space L2(M) and it is a legitimate question to ask what are their classical limits. Using 
the prescription ([M]) we obtain that: 

(P,{s,hr,,hO —-^h (37) 

which is trivially a function from Pi(R^). Hence, after performing the Fourier transform 
of ([371) . all the matrices |m)(m| are mapped in the classical limit to the same probability 
measure 6{q)6{p). Of course, one would expect that from the form of the energy spectrum 
of a harmonic oscillator, but as we said before, we consider the limit ([371) only as a formal 
exercise. The physically sensible classical limit of the Fock states is given by /i ^ 0, 
m — > oo, hm = const. In this limit one indeed recovers the classical microcanonical 
distribution function of the harmonic oscillator, as it was proven in Ref. [23] using the 
closely related method of Wigner functions. 

6 Concluding remarks 

The next logical step would be to try to apply the developed formalism to systems with 
compact kinematical groups, like, for example, spin systems with G = SU{2). The goal 
would be to describe the well known heuristic prescription: ^ — > 0, j ^ oo, jh = const, 
where j labels the irreducible representations of SU{2), within the presented group- 
theoretical formalism. One problem immediately arises: the corresponding classical phase- 
space is a sphere §^ [H] and Bochner's Theorem, crucial to our approach, holds only for 
Abelian groups. Thus, it is not so obvious what mechanism would allow one to recover 
classical statistics in this case. This is the subject of our further research (for an alterna- 
tive approach using co-adjoint orbits method see e.g. Ref. [24]; for another one based on 
non-commutative spheres see e.g. Ref. [25]). 

Another point is that at this stage our approach lacks a clear operational meaning of the 
mathematical concepts involved. Perhaps the most operationally flavoured reformulation 
of quantum theory is the one given by the quantum logic and orthomodular lattices (see 
Refs. [26l[27] for an introduction), as it operates directly with the probabilities of outcomes 
of (idealized) measurements. It also very coherently incorporates classical and quantum 
statistics within a common language. However, it comes with its own set of problems: the 
justification for the use of Hilbert spaces for building the lattice and the apparent lack of a 
clear quantum-to-classical transition mechanism within the formalism. Note that another 
problem: the justification for the lattice orthomodularity, was solved only very recently 
by Grinbaum [28] , using information-theoretical arguments. 

Summarizing, our work presents an alternative to the standard as well as to the alge- 
braic and lattice approaches to quantum statistics. It incorporates an elegant form of the 
quantum-to-classical transition. By the latter we mean a clear mechanism of showing how 
the classical state-space directly arises from the quantum one. 

We would like to thank V. Cappellini, J. Derezihski, J. Kijowski, A. Kossakowski, B. 
Mielnik, J. Naudts, K. Zyczkowski and especially M. Kus and M. Przanowski for dis- 
cussions, and the Deutsche Forschungsgemeinschaft (SFB 407, 436 POL), ESF PESO 
QUDEDIS, EU IP Programme "SCALA", and MEG (Spanish Goverment) under contract 
FIS2005-04627 and Gonsohder-Ingenio 2010 "QOIT" for the financial support. 
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